Friedmann cosmology with bulk viscosity: a concrete model for dark energy 
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The universe content is considered as a non-perfect fluid with bulk viscosity and can be described 
by a general equation of state (endowed some deviation from the conventionally assumed cosmic 
perfect fluid model). An explicitly bulk viscosity dark energy model is proposed to confront consis- 
tently with the current observational data sets by statistical analysis and is shown consistent with 
(not deviated away much from) the concordant A Cold Dark Matter (CDM) model by comparing 
the decelerating parameter. Also we compare our relatively simple viscosity dark energy model with 
a more complicated one by contrast with the concordant ACDM model and flnd our model improves 
for the viscosity dark energy model building. Finally we discuss the perspectives of dark energy 
probes for the coming years with observations. 

PACS numbers: 04.50.-|-h,04.25.Nx 



I. INTRODUCTION 

The cosmological observations indicate that the ex- 
pansion of our universe accelerates Recently lots 
of research work on its possible mechanism, such as ex- 
tended_gravity [2, decaying A Cold Dark Matter (CDM) 
model[3], modifying equation of state (hereafter EOS)or 
by introducing kinds of the so called dark energy mod- 
els are to explain the cosmic acceleration expansion ob- 
served. To consider causality and the hydrodynamical 
instability, an interesting barotropic dark energy model 
is proposed^ that includes a linear EOS of with a gen- 
eral form, p — a{p — pq) which incorporated into 
cosmological model can describe the hydrodynamically 
stable dark energy behaviors. 

The astrophysical observations also indicate that the 
universe media is not a perfect fluid [6,] and the viscosity 
is concerned in the evolution of the universe 0, B Hi ■ On 
the other hand, in the standard cosmological model, if 
the EOS parameter uj is less than —1, the universe shows 
the future finite singularity called Big Rip [l^, [HI • Sev- 
eral ideas are proposed to prevent the big rip singularity 
as thought it un-physical, like by introducing quantum 
effects terms in the action 12], or by including universe 
viscosity media for the Universe evolution [13]. So an in- 
teresting question naturally arises: what kind of role the 
cosmic viscosity element can play for helping the above 
two facets, dark energy and cosmic dark energy model 
singularity? Considering some deviation from the ideal 
fluid model is also very helpful to nowadays cosmology 
probes advancement. 

In Refs. [il [H, [H, [O, d, d, HO], the bulk viscosity 
in cosmology has been studied in various aspects. Dissi- 
pative processes are thought to be present in any realistic 
theory of the evolution of the universe. In the early uni- 
verse, the thermodynamics is far from equilibrium, the 



viscosity should be concerned in the studies of the early 
stage for cosmological evolution. But even in the later 
cosmic evolution stage, for example, the temperature for 
the intergalactic medium (IGM), baryonic gas, generally 
is about lOOOOK to lOOOOOOK and the complicated IGM 
is rather non-trivial. The sound speed Cg in the bary- 
onic gas is only a few kms~^ to a few tens kms~^ and 
the Jeans length A yields a term as an effective viscosity 
CsX. On the other hand , the bulk velocity of the bary- 
onic gas is of the order of hundreds kms~^ [2l|. So it is 
helpful to consider the viscosity element in the later cos- 
mic evolution. It is well known that in the framework of 
Friedmann Robertson Walker (FRW) metric, the shear 
viscosity has no contribution in the energy momentum 
tensor, and the bulk viscosity behaves like an effective 
pressure. At the late times, since we do not know the na- 
ture of the universe content (dark matter and dark energy 
components) very clearly, concerning the bulk viscosity 
is reasonable and practical. Moreover, the cosmic viscos- 
ity here can also be regarded as an effective quantity as 
caused by complicated astrophysics mechanisms and may 
play a role as a dark energy candidate [12, IH, 1131 ^ 
possible unification scheme for the two mysterious dark 
components (dark matter and dark energy) [U [H, [13] 
as they may be facets of the same problem, and then 
evidence for dark matter is also evidence for dark energy. 

The EOS p — — l)p and the bulk viscosity ^ — ap'^ 
is studied in the full causal theory of bulk viscosity, and 
the case s = 1/2 has possessed exact solutions as shown 
in [Uj [13] ■ However, both the pressure and the bulk 
viscosity coefficient may have constant components. We 
argue that the non-causal approximation is reasonable 
in the late times of the universe evolution. In our previ- 
ous papers, we have shown that the Friedmann equations 
can be solved with both a more general EOS and bulk 
viscosity detailed as follows. (Note that the EOS with a 
constant pressure gets degenerate with ACDM model.) 
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P=(7-l)p + Po, (1) 
where po and 7 are two parameters. The bulk viscosity 
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is expressed as 



C = Co + Ci-- 

a 



(2) 



where Co and C,i are two constants conventionally, and 
the overhead dot stands for derivative with respect to 
time. The motivation of considering this bulk viscosity is 
that by fluid mechanics we know the transport/viscosity 
phenomenon is involved with the "velocity" a, which is 
related to the scalar expansion 6 = 3d/ a. Both ^ = 
Co (constaiit) and Q oz 6 are considered in the previous 
papers 0, |2g|, so a linear combination of the two are 
more general. The to = p/ p is constrained as —1.38 < 
uj < —0.82 by present observation data, so the inequality 
in our previous case should be 



- 1.38 < 7 - 1 + — < -0.82. 
P 



(3) 



The parameter po can be positive (attractive force) or 
negative (repulsive force), and conventionally Co and Ci 
are regarded as positive. To choose the parameters prop- 
erly, it can prevent the Big Rip problem or some kind of 
singularity for the cosmology model, like in the phantom 
energy phase, as demonstrated before. Additionally, the 
sound speed in this model can also keep the casuality 
condition. 

This present paper is a continuous work following our 
previous efforts, which is organized as follows. In Sec. 
II we present a relatively simple cosmology model with 
extremely non-relativistic dark matter and viscosity dark 
energy by an explicit form. With these we give out the 
exact solution and discuss the acceleration phase in this 
model. In the last section (Sec. Ill ) we discuss and sum- 
marize our conclusions. 



II. MODEL AND CALCULATIONS 

We consider the Friedmann-Roberson- Walker metric 
in the flat space geometry (A;=0) as the case favored by 
WMAP satellite mission on cosmic background radiation 
(CMB) data 



ds^ = -dt^ + a{tf{dr^ + r^dn^), 



(4) 



and assume that the cosmic fluid possesses a bulk viscos- 
ity. The energy-momentum tensor can be written as 



T^u = pUf,U, + ip + n)H, 



(5) 



where in the co-moving coordinates L/^ — (1,0), and 
Hfj,i, = g^ii, + UfiUv [25i]. By defining the effective pres- 
sure as p = p -\- Ti and from the Einstein equation 
-R/^f ~ ^gfiiyR ~ 87rG'r^i^, we obtain the Friedmann equa- 
tions 

d2 8ttG 



AttG 



{p + 3p). 



(6a) 
(6b) 



The covariant conservation equation for energy T^J^, 
yields 



p+ip + p)9 ^ 0, 

where the expansion parameter 9 — Ufjj^ = 3d /a. 
Then the Friedmann equations give: 



(7) 



SttGp = 3 



= 3H^ 



If we assume that the cosmic fluid possesses a bulk vis- 
cosity as shown explicitly in the pressure as 11 = —(^9 
(we take another Greek character to represent the bulk 
viscosity which is different from the before papers), the 
energy-momentum tensor could be written fully as : 

T^. - pU^U, + (p - i9) H^, 

where denotes 4 - velocity, as before defined 9 = 3H, 
and ^ is the bulk viscosity with its explicit form to be 
present below. 

The cold dark matter has been assumed to be ex- 
tremely non-relativistic so that we can take its pressure 
p = 0, and we now suppose that the effect of dark en- 
ergy on cosmos evolution is included in the viscous term 
—^9 which has the dimension of pressure. In this present 
work we treat it as the effective pressure of dark energy. 
Then the covariant conservation of T^j, (Eq.(7)) yields: 



p+ip-^9)9^0 



so, 



P 



= S_9^ 



(8) 



Here, it will be convenient to define a dimcnsionlcss pa- 
rameter below: 

Hq Per 

where pcr is the critical density. Using this definition, we 
could rewrite equation ([5]) as: 



dh'^/dt 
Ho 



3h'^ = 9Xh^ 



(9) 



where A = £,Ho/pcr is the bulk viscosity. 

For comparison with the observational data, we should 
derive h from Eq. ^ by writing it as the equation of red- 
shift z as: 



2(1 + z)^ + 3h = 9\ 

dz 



(10) 



In the last part of [33| , condition with constant viscosity 
has been considered. It proves that it has a good fitting 
with SNe la data sets. Recently, a complicated viscosity 
form dependent on h has been proposed in ref.js^. To 
overcome old cosmological constant problem we believe 
that a variable viscosity will give a better result. Here 
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we propose an explicit red-shift dependent bulk viscos- 
ity, which gets a constant limit (an effective cosmological 
constant) today as z = 0. And we will show that this 
relatively simple form is better than the complex one in 
ref. [3^] when compared with ACDM model. The explicit 
red-shift dependent viscosity is 



9A = Ao + Ai(l + z)'^ 



(11) 



where n is an arbitrary integer, Aq, and Ai arc two arbi- 
trary constants, which could all be best fitted from the 
observational data sets. 

With this new bulk viscosity, we could write Eq. pOjl 
as: 

-2(l + z)^ + 3/i = Ao + Ai(l + z)" (12) 
az 

which is a first-order differential equation of h and the 
exact solution of h from this equation is: 

/^ = A.(l + .)--^(l + z)" + ^ (13) 

where A2 is an integration constant. The obtained rela- 
tion above could be regarded that the cosmic expansion 
rate is from a combined result of the viscosity term: the 
matter component effect plus the effective dark energy 
(including the cosmological constant like term ^) con- 
tributions. Because of the consistent requirement h = 1 
for the spatial fiat universe at present, from the above we 
have 

3 2n - 3 

With the analysis of dimension of term A2 and the fact 
that matter density of our cosmos Pm oc a{t)^^, where 
a{t) is the scalar scale factor we will show that the value 
of A2 from data fitting below does consist with the result 
of matter component given by WMAP5 data setsjll]. To 
make the fitting result compared with ACDM model, we 
make an assumption that A2 = ^/flmo, where ftmo = 
%ttG pmo / {^Hq^) . So in our viscosity dark energy model, 
the first term as labelled by A2 could be looked as mat- 
ter contribution and the other terms of h{z) is due to 
the viscosity dark energy, which with limit case {z — Q 
today) as the concordant ACDM model. The consistent 
requirement relation also gets clear meaning, that is, the 
relationship between the matter and dark energy (includ- 
ing the cosmological constant contribution) components. 

The observations of the SNe la have provided the first 
direct evidence of the accelerating expansion for our cur- 
rent universe, so any model attempting to explain the ac- 
celeration mechanism should be consistent with the SNe 
la data implying results, as a basic requirement. As we 
know the observations of supernovas measure essentially 
the apparent magnitude m, which is related to the lumi- 
nosity distance by 



m 



M + 51ogio Dl [z) 



(14) 



where the distance Dl (z) = {Hq) (z) is the dimen- 
sionless luminosity and 



= (1 + z)dM{z), 
where d^/ is the co-moving distance given by 



1 



H{z') 



dz 



Also, 



M = M 



5 logio I I 



\lMpc J 



25, 



(15) 



(16) 



(17) 



where M is the absolute magnitude which is believed 
to be constant for all supernovaes of type la. In this 
paper, we use the 307 Union SNe la data sets compiled 
in reference [26] . The data points in these samples are 
given in terms of the distance modulus 



p,obs = m{z) - Mobs{z) 



(18) 



We employ it for doing the standard statistic analysis. 
So the is calculated from 



5 logio D 



Lth 



{Zi, C^^ ) 



<7obs{zi) 



(19) 



where M — Ai — Mobs is a free parameter and 
DLth{zi]Ca) is the theoretical prediction for the dimcn- 
sionless luminosity distance of a supernovae at a partic- 
ular distance, for a given model with parameter Ca- 

On the other hand, the shift parameter TZ and the dis- 
tance parameter A are considered to give contributions to 
data fitting. The shift parameter TZ is defined in refs. [2§| 
and 1301 as 



n 




dz 

W) 



(20) 



and WMAP5 results [281 have updated the redshift of 
recombination to be z^: = 1090. Its detail meaning can 
be found in reference [3l|. The distance parameter A is 
defined as 



A: 



HzbY 



, 1 

' Zb Jo 



dz 

W) 



(21) 



where Zb — 0.35. 

Considering TZ and A, we use the total to make the 
standard statistic analysis, and data fitting: 



^total 



X 



n-TZob 

o■7^ 



(22) 



The resuh of data fitting is fisted in TABLE 1. 
In Fig. 2, we plot the deceleration parameter q relation 
with the redshift as model comparisons: 



1 



H 



z dH 
dz 



1 



(23) 
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FIG. 1: The la, 2a, 3a C.L. contours of Ai and A2 in the best 
fitting condition n — —1. The black dot corresponds to the 
best fitting values. 



TABLE I; Fitting results for model parameters 



n 


Ao 


Ai 


A2 


x' 


-1 


0.18027 


2.0923 


0.52145 


312.8 


-0.8 


-0.04672 


2.2784 


0.52027 


313.0747 


-2 


0.63824 


1.8331 


0.52538 


312.2193 


1 


1.0464 


-0.15136 


0.50016 


323.3296 



At the same time, we plot the deceleration parameter 
relations of q{z) in ACDM model and the more complex 
bulk viscosity in ref.[^] for contrast in Fig. 2, too. 

To demonstrate clearly how and when the viscosity 
dark energy components catch up of the matter con- 
tribution, and overpass it we also plot the two compo- 
nents evolution in Fig. 3 With our best fitting result of 
the n = —1 case, we could calculate the EOS of viscosity 
dark energy. The pressure of dark energy is —^0, so state 




FIG. 2: q — z relation. Red line (the lowest one) represents 
ACDM, blue thick line and black line are our models with n — 
— 1 and n = 1 respectively. By including the more complex 
viscosity form in [32], we get another result (dashed line) for 
illustration. 



FIG. 3: The bulk viscosity dark energy part (the slow chang 
ing almost horizontal line) evolutions vs matter component 



parameter: 

PX 

where px = p(l ~ f^m) is the density of dark energy. 
Then we have, 

3iJ^ 



(24) 



p{l - 0„) 



(25) 



With the relation of ^ and A defined above, we could get 
Lo today: 



cjo = - 



Aq + Ai 

3(1 — ilmo) 



(26) 



By the best fitting data we give an approximated value 
of EOS at present(/i = 1) cj ~ -1.04. There is a little 
deviation from ACDM model by a value of -0.04. 

Recently, ref.jl^l has proposed Ora diagnose method 
to differentiate a new model from the ACDM model with 
the constant equation of state parameter exactly as -1. 
The diagnostic parameter Om is defined as: 



Cm = 



(27) 



where x = \ -\- z. The dimensionless expansion param- 
eter h(z) for a dark energy model with constant equation 
of state can be written clearly as: 



h{xf = f^„o^' + (1 - f^™o)2:'^'+"^ 



(28) 



So, for the ACDM model with its EOS parameter uj — 
— 1, Om = Omo- This result also gives us a null — test of 
cosmology constant. For model we proposed(here we use 
the best fitting results with n = — 1, which gives a better 
q{z) evolution compared with the ACDM model), its 
could be written as: 



= (A: 



Ai 



2X ' + —X 



3 ' 



(29) 



so the Om evolution for this viscosity dark energy model 
is shown in Fig. 4 to compare with the ACDM model 
which is a horizontal line in the plot. 
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FIG. 4: Om — z relation as clearly shown the viscosity dark 
energy model deviations from ACDM model. The dashed line 
corresponds ACDM model while the thick curve corresponds 
to our viscosity dark energy model contribution. 
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FIG. 5; q{z) — 2 relations for three different models. The 
nearby two lines are ours with the explicit viscosity form and 
the concordant ACDM model(f2 — 0.3). The dashed line 
corresponds the complicated one, which can degenerate the 
previous two models for certain redshift values 

With this diagnostic method, we could also see that 
this model has different properties from ACDM model 
with the exact constant EOS parameter lu = — I. 

With the above discussions we can see our relative sim- 
ple viscosity model fits the concordant ACDM model bet- 
ter. It is due to the simpler viscosity form. As a matter 
of fact so far we still have not obtained a systematic way 



to express the cosmic viscosity effects. In view of the 
astrophysical reality any progresses in this field is very 
helpful. To show the similarities and differences among 
the concordant ACDM model, our model and the com- 
plicated one [3^ we present a much clear figure below 
with the hope there will be much relative works to ap- 
pear soon. Hope the study for viscosity dark energy can 
be stimulated further. 



III. SUMMARY AND DISCUSSIONS 

In this letter we present an explicit viscosity form to 
mimic dark energy behaviors and confront it with cur- 
rent observational data sets. Though the viscosity dark 
energy form is simple, the results are better. The possi- 
ble Universe evolution fate, like possible future singular- 
ity types from this viscosity model can be also discussed 
and we will work it out elsewhere. We emphasize that 
perfect fluid is just a limit case of a general viscosity me- 
dia that is more practical in the astrophysics sense. 

Discovery of dark energy is about ten years old, but its 
nature and origin have been still puzzling. Fundamental 
as it has promised to physics foundations there are sev- 
eral possibilities to develop with great expectations. One 
of these is the unification scenario of dark matter and 
dark energy, that is, they are two facets of one secret. If 
we finally discover the dark matter either by accelerators 
such as LHC or future ILC or by satellite missions, like 
PAMELA and GLAST now on the running, we may in- 
fer the existence of dark energy therefore. In this aspect 
viscosity dark energy has already shown its unification 
efforts, to write uniformly all the possible cosmic compo- 
nents in a formula with distinct scaling law with respect 
to scalar factor [ssl. [s^ . It is worthy further studying. 
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